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Abstract  -  In  target  positioning  and  tracking,  most 
sensors  provide  measurements  either  as  range  or  bearing 
or  both.  The  measurements  are  used  to  update  an  a  priori 
estimate  either  via  a  linearized  least  squares  method  or  an 
extended  Kalman  filter.  In  either  case,  the  resulting 
solution  has  two  components,  one  is  related  to  the 
measurement  prediction  errors  and  the  other  is  an 
observation  matrix  obtained  from  linearizing  the 
nonlinear  measurement  equations  around  the  a  priori 
estimate.  This  paper  studies  the  geometric  factors 
explicitly  and  relates  the  observation  matrix  to  the  line  of 
sight  (LOS)  vector  for  a  ranging  sensor  and  the  direction 
perpendicular  to  the  LOS  vector  of  a  bearing-only  sensor. 
As  a  result,  the  updating  of  estimation  error  covariance 
with  range  and  bearing  measurements  can  be  intuitively 
assessed  via  the  shaping  of  estimation  error  ellipse  along 
LOS  directions.  It  provides  a  valuable  means  for  target 
positioning  and  tracking  performance  modeling  and 
prediction  and  can  thus  be  used  in  active  management  of 
distributed  sensor  resources  and  sensor  path  planning. 
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1  Introduction 

Passive  and  active  sensors  are  widely  used  in  target 
position  location  and  tracking  [1,  2,  3,  7,  10].  Passive 
sensors  include  (1)  acoustic  sensors  that  detect  ground 
vehicle  motion  and  vibration  and  (2)  imaging  sensors  that 
capture  reflected  visible  lights  or  radiated  thermal  energy 
from  a  target.  In  addition  to  signatures  characteristic  to 
targets  for  possible  identification,  these  passive  sensors 
also  measure  the  direction  of  arrival  of  the  mechanic  or 
electromagnetic  signals  of  the  targets  in  terms  of  their 
bearing  angles  relative  to  the  sensor  platforms  for 
positioning  and  tracking.  Additional  examples  include 
passive  sonar  sensors  for  underwater  applications  and 
passive  radar  sensors  for  space  observations  [8]. 

Active  sensors  include  radar,  laser  radar,  and  active  sonar. 
By  measuring  the  round  trip  time  of  flight  of  an  energy 
pidse,  active  sensors  provide  range  measurements  to  targets. 
Time  difference  of  arrival  (TDOA)  is  a  common  technique 
for  differential  range  calculations.  Some  active  sensors  also 
measure  the  angles  of  arrival  (AOA)  of  returned  pulse  or 
the  direction  in  which  the  transmit  beam  is  pointed. 
Monopulse  and  antenna  array  techniques  provide  more 
accurate  angular  measurements  than  mechanically  scanned 
antennas.  Other  active  sensors  are  capable  of  extracting  the 
Doppler  frequency  shifts  impacted  by  target  motion,  thus 
measuring  the  range  rate. 


According  to  the  type  of  measurements,  most  sensors  used 
in  target  positioning  and  tracking  can  be  viewed  as 
measuring  either  range  or  bearing  or  both  (without 
considering  range  rate  and  target  signature  measurements) 
[10].  A  target’s  position  is  determined  using  multilateration 
with  ranging  sensors  or  triangulation  with  angular  sensors 
or  their  mixture.  Multilateration  or  triangulation  can  be 
done  when  a  target  is  viewed  by  multiple  distributed 
sensors  or  by  a  single  sensor  moving  along  (multiple 
views).  A  batch  processing  mode  may  use  the  conventional 
least  squares  (LS)  method  whereas  a  sequential  processing 
mode  calls  for  the  Kalman  filter  [1,  2,  3,  7,  10]. 

Indeed,  the  sensor  measurements  are  used  to  update  an  a 
priori  estimate  either  via  a  linearized  least  squares  method 
or  an  extended  Kalman  filter.  In  either  case,  the  resulting 
solution  has  two  components,  one  is  related  to  the 
measurement  prediction  errors  and  the  other  is  an 
observation  matrix  obtained  from  linearizing  the  nonlinear 
measurement  equations  around  the  a  priori  estimate.  The 
observation  matrix,  which  determines  the  estimation  error 
covariance,  has  a  distinct  structure  containing  the 
geometric  information  for  a  given  target  positioning  and 
tracking  setting. 

It  is  well  known  that  the  observation  matrix  is  related  to  the 
sensing  geometry.  However,  this  paper  elucidates  the 
geometric  factors  for  both  active  and  passive  sensors  [4,  14, 
15]  and  explicitly  relates  the  observation  matrix  to  the  line 
of  sight  (LOS)  vector  for  a  ranging  sensor  and  the  direction 
perpendicular  to  the  LOS  vector  of  a  bearing-only  sensor. 
The  geometric  expression  provides  a  mathematical  basis 
for  the  well-known  fact  that  the  error  ellipses  for  co-located 
passive  and  ranging  sensors  are  perpendicular. 

Furthermore,  the  updating  of  estimation  error  covariance 
(i.e.,  the  eigen  structure)  with  range  and  bearing 
measurements  can  be  intuitively  assessed  via  the  shaping  of 
estimation  error  ellipse  along  LOS  directions.  For  an 
anticipated  tracking  scenario,  this  provides  an  efficient 
means  for  target  positioning  and  tracking  performance 
modeling  and  prediction  and  can  thus  be  used  in  active 
management  of  distributed  sensor  resources  and  sensor 
path  planning  [5,  6,  9,  16]. 

The  rest  of  the  paper  is  organized  as  follows.  In  Section  2, 
the  least  squares  method  is  used  to  relate  the  position 
solution  to  the  observation  matrix  and  to  the  LOS  vectors. 
In  Section  3,  the  updating  of  error  covariance  by  LOS 
vectors  is  characterized  with  simulation  examples.  Finally, 
the  paper  is  concluded  in  Section  4. 
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2  Geometric  Scaling  via  LOS  Vectors 

In  this  section,  we  first  present  the  position  solution  via  the 
least  squares  method.  We  then  relate  it  to  the  observation 
matrix  via  linearization.  Next  we  express  the  observation 
matrix  in  terms  of  LOS  vector  for  active  and  passive 
sensors. 

2.1  Least  Squares  Solutions 

Consider  a  target  at  x  and  the  i-th  sensor  at  x,-.  The  i-th 
sensor’s  measurement  is  given  by: 

A  =  /(*,*,)+v,  (!) 

where  /X-,  •)  is  a  nonlinear  measurement  equation  and  v,  is 
the  sensor  measurement  error  being  zero-mean  Gaussian 

Ho,  <r2). 


x  =  argmin(z  -  Hx)7R  '(z  -  Hx)  +  (x-x0)7P0  ‘(x  -x0)  (7) 

X 

The  least  squares  solution  is  given  by: 

x  =  (HrR 'H  +  P0-1)-1(HrR"‘z  +  P0“‘x0)  (8a) 

P  =  (HrR  H  +  P0)  1  (8b) 

Alternatively,  we  can  treat  the  a  priori  knowledge  x0  and 
P0  as  an  additional  measurement  and  obtain  the  augmented 
measurement  equations  as: 

z  =  Hx  +  v  (9a) 

x0=x  +  £  (9b) 

v  ~  HO,  R),  e  ~  HO,  P0),  E  { vs7’}  =  0  (9c) 

Now  apply  the  Kalman  filter’s  measurement  updating  step 
to  the  augmented  measurement  [zr  x„]r  as: 


Denote  an  initial  estimate  of  the  target  state  by  x0.  The 
nonlinear  measurement,  be  a  range  or  a  bearing,  can  be 
linearized  around  the  estimate: 
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In  terms  of  measurement  prediction  error,  the  equation  can 


be  further  written  as: 

z,  =z;-/}(Xo,x,.)  +  hfxo  =hfx  +  v;  (3) 

For  m  sensors,  the  linearized  measurements  can  be  put  into 
a  vector  format  as: 

z  =  Hx  +  v  (4a) 

2T  =[zt  z2  •••  zj  (4b) 

vr=[v,  v2  •••  vj  (4c) 

Hr  =  [h,  •••  hj  (4d) 

An  estimate  of  the  target  state  given  the  measurements  can 
be  obtained  as  the  least  squares  (LS)  solution  to  the 
following  performance  index  as: 

x  =  argmin(z-Hx)rR_1(z-Hx)  (5a) 

x 

where  R  is  the  measurement  noise  covariance  matrix 
defined  as: 

R  =  £{vv7}  (5b) 

The  least  square  solution  is  given  by: 

x  =  (HR  H)  HR  z  (6a) 

P  =  £{(x-x)r(x-x)}  =  (H7R  ’H)"1  (6b) 


Assume  that,  in  addition  to  the  initial  estimate  of  the  target 
state  x0,  the  estimation  error  covariance,  denoted  by  P0,  is 
also  known.  Then,  given  the  target  measurements  as  in  (4a), 
we  can  obtain  the  weighted  least  square  (WLS)  solution 
that  minimizes  the  following  performance  index  as: 


The  terms  in  front  of  [zr  x7]7  in  (10)  constitute  the 

Kalman  filter  gain.  Expanding  the  matrices  in  terms  of  their 
elements  leads  to  the  solution  (8a).  Similarly,  the  Kalman 
filter’s  covariance  updating  step  can  be  written  up  for  the 
augmented  equations,  leading  to  the  updated  covariance 
(8b).  Note  that  the  inverse  of  the  covariance  matrix  in  (8b) 
is  also  called  the  information  matrix,  denoted  by  A: 

A  =  P  1  =HrR  ‘H  +  P0  ‘  (11) 


2.2  Active  Ranging  Sensors 

First  consider  the  case  with  active  sensors  (ranging 
measurements)  in  a  two  dimensional  setting.  The  nonlinear 
range  equation  (1)  is  spelled  out  as  (ignore  the  noise  term): 

r,  =  f ,{*,*,)  =  V(*  -  x,)2  +  (y  -  y,-)2  (12a) 


The  corresponding  linearized  equation  around  an  initial 
estimate  x0  of  x  is: 

r,  “/•(x0,xi)  +  hf(x-x0)  (12b) 


with 
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It  is  clear  from  (12b)  that  h;  is  the  line  of  sight  (LOS) 
vector  from  the  ;-th  sensor  to  the  target,  denoted  by  e;,  as 
shown  in  Fig.  1. 


2.3  Passive  Bearing-Only  Sensors 

Referring  to  Fig.  1,  we  now  consider  the  case  with  passive 
sensors  (bearing  measurements)  still  in  a  two  dimensional 
setting. 
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The  nonlinear  angular  equation  is: 

/,(*,*,)  =  9,  =  tan-1  (— — — )  (13a) 

X  —  X- 

which  can  be  linearized  around  an  initial  estimate  x0  of  x 
as: 


3  */Xxo>xi)  +  hf(x-xo) 


(13b) 
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Comparing  (13c)  to  (12c),  it  shows  that  the  measurement 
matrix  h,  for  the  bearing-only  sensor  is  range-dependent. 
More  importantly,  it  is  perpendicular  to  the  line  of  sight 
vector  from  sensor  to  target.  When  the  range-dependence 
and  angular  errors  are  combined,  it  provides  a  position 
error  of  r,cr,  along  the  direction  perpendicular  to  the  LOS 
(i.e.,  along  e,1). 


2.4  Geometric  Dilution  of  Precision  (GDOP) 

A  scalar  value  that  characterizes  the  position  solution  is  the 
geometrical  dilution  of  precision  (GDOP)  defined  as  [4, 
13]: 

GDOP  =  •y/trace((HrR_1H)“1)  (14a) 

=  A/trace((HrH)-‘)  when  R  =  I  (14b) 

For  the  case  with  two  ranging  sensors,  the  GDOP  can  be 
written  as: 
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For  the  case  with  two  bearing-only  sensors,  the  GDOP  can 
be  written  as: 
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Eqs.  (15)  and  (16)  are  consistent  with  our  intuition  that  two 
perpendicular  sensors  can  produce  the  best  estimate  and 
two  collinear  sensors  cannot  produce  a  solution.  Different 
from  ranging  sensors,  however,  bearing-only  sensors  have 
range  terms  explicitly  as  in  (13c).  In  fact,  the  error  ellipses 
for  passive  and  ranging  sensors  are  perpendicular  are  well 
known.  It  is  meaningful  to  explicitly  show  it  from  the 
observation  matrix  in  terms  of  LOS  vectors. 

As  shown  next,  the  line  of  sight  vectors  and  measurement 
error  covariance  matrices  are  related  to  the  eignevalues  and 
eigenvector  directions  of  the  estimation  error  ellipses. 
When  used  to  update  a  position  solution,  they  shape  the 
error  ellipses  by  adjusting  its  size  and  orientation  via  the 
eigenstructure. 

3  LOS  Updating  and  Characterization 

As  shown  in  the  last  section,  the  linearized  measurement 
equation  for  an  active  ranging  sensor  or  a  passive  bearing- 
only  sensor  is  characterized  by  (1)  the  line  of  sight  (LOS) 
vector  or  its  perpendicular  and  (2)  the  measurement  error 
variance.  A  position  solution  is  also  determined  by  these 
two  factors  as  is  evident  from  the  following  equation  of 
estimation  error  covariance: 

P  =  (HrR  'H  +  Po"1)"1  (17) 

where  P0  is  the  a  priori  covariance  matrix,  P  is  the  a 
posteriori  covariance  matrix,  H  is  the  measurement  matrix, 
and  R  is  the  measurement  error  covariance. 

To  simplify  the  analysis,  first  consider  the  case  where  Po  = 
I,  H  =  ar  =  [cos#  sin#],  which  is  a  row  vector  with  #being 
the  angle  relative  to  the  x-axis,  and  R  =  cf .  This  represents 
the  case  where  the  a  priori  estimate  is  updated  by  a  single 
measurement.  The  updated  covariance  is  given  by: 

(HrR  'H  +  P;1)"1  =(I  +  o-"2aar)“‘  =  l-yaaT  (18a) 


where 

r_  1  _  1  -f1’  a  =  0  (18b) 

a'a  +  cr2  1  +  cr2  [o,  <t  =  co 

It  is  easy  to  verify  that  a  is  an  eigenvector  of  the  matrix  (I  - 
raa T)  with  the  corresponding  eigenvalue  (1  -  yaTa)  =  1  -  y. 
For  the  2D  case,  the  other  eigenvector  is  perpendicular  to  a, 
denoted  by  a1,  with  the  corresponding  eigenvalue  being 
unity. 
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In  terms  of  geometry  interpretation,  the  unity  circle  is 
shaped  into  an  ellipse  with  the  following  eigenvalues  and 
eigenvectors: 


A  =  diag([l-y  1]) 
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The  following  simulation  is  used  to  illustrate  the 
characteristics  of  LOS  updating.  Let  the  initial  covariance 
be  a  unit  matrix,  P0  =  diag([  1,  1]),  which  represents  a 
circular  error  as  shown  in  Fig.  2.  We  vary  the  LOS  vector’s 
orientation  6  relative  to  the  x-axis  and  the  measurement 
error  covariance  R  o2.  Fig.  2  shows  two  updates  at  (#,  a) 
=  (30°,  1)  and  (70°,  0.5),  respectively.  Each  LOS  update 
squeezes  the  circle  into  an  ellipse  with  different  orientation 
and  different  eccentricity.  The  orientation  is  such  that  the 
eigenvector  corresponding  to  the  small  eigenvalue  points 
along  the  LOS  direction. 

Fig.  3  shows  eigenvector  constriction  where  the  major  and 
minor  eigenvectors  are  off  by  90°  as  expected  and  the 
minor  eigenvector’s  angle  (j)  changes  linearly  as  the  LOS 
angle  #. 

The  major  eigenvalue  remains  constant  but  the  minor 
eigenvalue  changes  inversely  with  a  (according  to  (15b)) 
as  shown  in  Fig.  4.  The  better  the  sensor  performance  (the 
smaller  a),  the  smaller  is  the  minor  eigenvalue  after  update 
and  vice  versa. 

If  a  sensor  stares  at  a  target  in  the  same  direction,  the  error 
ellipse  is  updated  repeatedly  by  the  same  LOS  vector 
(assume  the  target  is  not  moving).  Then  the  reduction  of  the 
minor  eigenvalue  over  time  (i.e.,  the  number  of  updates)  is 
shown  in  Fig.  5  for  different  sensor  quality  (i.e.,  R).  This  is 
equivalent  to  time  averaging. 

For  a  single  sensor,  another  way  to  reduce  the  errors  in 
target  positioning  is  to  circle  the  target  as  in  a  spot 
synthetic  aperture  radar  (SAR)  operation.  Update  in  one 
direction  squeezes  the  error  in  that  direction.  When  it  is 
done  in  all  directions,  it  makes  a  smaller  error  circle  as 
shown  Fig.  6  where  the  errors  at  5  discrete  angles  are 
shown  together  with  the  initial  large  circle. 

Fig.  7  shows  the  reduction  of  the  eigenvalues  as  a  function 
of  LOS  angle  0.  As  shown,  the  major  eigenvalue  decreases 
rather  quickly  over  the  first  90°  and  then  much  slowly  for 
next  90°  and  does  not  change  much  for  the  remaining  180°. 
This  indicates  that  a  single  ranging  sensor  can  significantly 
reduce  positioning  error  by  circling  the  target  by  90°  or  so, 
which  is  consistent  with  our  intuition  and  can  be  used  as  a 
practical  guideline  in  sensor  placement  and  scheduling. 

Consider  the  case  with  two  LOS  updates  simultaneously. 
The  two  sensors  have  their  LOS  vectors  sl\T  =  [cos#i,  sin#i] 
and  a2r  =  [cos#;,  sin#>]  with  erf  and  os2,  respectively.  The 
updated  covariance  is  given  by  [8]: 


(I  +  oya.af  +cr:2a,a()“1  = 

=  I-Kw»f  +72*2*1  -PTxYi^A -PYiYi*2*[)  (20a) 

where 

P  =  afa2  (20b) 

v,= _ ! _ =  (20c) 

of  +  afa,  of+l 

v,= _ I _ =  _ L_  (20d) 

/  2  2  T  2  i 

O’ 2  "^”^2^2  £7 ’2  “I-! 

r  = _ I _  (20e) 

l-fViiU 


Now  consider  two  special  cases.  In  the  first  special  case,  a; 
=  a2  but  there  is  unequal  variance  erf  ^  os2.  That  is,  a  target 
is  detected  in  the  same  direction  but  at  different  ranges.  It  is 
easy  to  verify  that  (20a)  becomes: 


(I  +  a-Vfr^I-^af 

(21a) 

where 

a "2  =  er,“2  +  ay 

(21b) 

_  _  i  cr 2  +  a\ 

Y  a2  +1  of  +  of  +  of  cr2 

(21c) 

If  erf  =  <r22,  then  a2  =  of !  2  for 

an  additional  update.  This 

is  equivalent  to  reducing  the  measurement  error  variance 
by  n  with  n  additional  updates.  The  corresponding  eigen 
value  decreases  according  to  l/(/7+l)  for  unity  variance  (erf 
=  1).  The  evaluation  of  continuously  updated  covariance  is 
similar  to  running  the  covariance  update  equation  of  a 
Kalman  filter  (i.e.,  the  Riccatti  equation). 

In  the  second  special  case,  a!  _L  a2,  that  is,  they  are 
perpendicular  to  each  other.  Then  p  =  0  in  (20b)  and  y  =  1 
in  (20e)  and  (20a)  becomes: 

(I  +  oya,af  +  aya,a()_1  =I-/,a,af  -f2a2a2  (22) 


It  is  easy  to  verify  that  (22)  has  the  following  eigenvalues 
and  eigenvectors: 

A  =  diag([\  -  yt  1  -y2J)  (23a) 

V  =  [a,  a2]  (23b) 

This  indicates  that  for  two  orthogonal  LOS  vectors, 
simultaneous  update  with  two  orthogonal  sensors  is 
equivalent  to  two  independent  updates  in  terms  of  the  final 
error  ellipse  and  their  eigen- structure. 

For  the  general  case  with  two  simultaneous  LOS  updates, 
the  updated  covariance  is  given  by  (20).  The  calculation  of 
its  eigenvalues  and  eigenvectors  is  rather  involved  except 
for  the  two  special  cases  given  in  (21)  and  (22).  The 
following  simulation  is  used  to  analyze  the  general 
eigenvalues  and  eigenvectors. 

Again  let  the  initial  covariance  be  a  unit  matrix,  P0  = 
diag([  1,  1]).  We  fix  the  first  LOS  vector  at  9\  =  30°  with  Rj 
=  1.  We  then  vary  the  second  LOS  vector  #2  from  0°  to 
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Fig.  2  Update  of  Unit  Circle  with  Different  cr  and  6 


Fig.  3  Eigenvector  Directions  ^vs.  LOS  Angle  6 


error  ellipse  re-sized  by  los  update 


update  quality,  a2 


error  ellipse  re-sized  by  repeated  updates 


Fig.  4  Eigenvalues  vs.  Sensor  Quality  <j 


Fig.  5  Eigenvalues  vs.  Number  of  Updates 


error  ellipse  re-shaped  by  continuous  circular  update 


error  ellipse  re-shaped  by  continuous  circular  los  update 


Fig.  6  Errors  over  Circular  Updates 


Fig.  7  Eigenvalues 
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180°  and  R2  from  0.1  to  10.  The  major  and  minor 
eigenvalues  as  a  function  of  ft  for  different  R2  are  showed 
in  Figs.  8  and  9,  respectively.  Since  ft  =  30°,  the  max  and 
min  values  for  both  major  and  minor  eigenvalues  appear  at 
02  =  30°  and  62  =  120°,  respectively,  which  represent  the 
two  special  cases  given  in  (18)  for  A9  =  62  -  9\  =  0°  and  in 
(19)forA6>=ft-ft  =  90°. 

Again  with  the  first  LOS  vector  fixed  at  ft  =  30°  with  Ri  = 
1,  Figs.  10  and  11  show  the  major  and  minor  eigenvector 
angles  as  the  second  LOS  vector  ft  varies  from  0°  to  180° 
and  R2  from  0.1  to  10.  The  two  eigenvector  angles  have  the 
same  shapes  except  that  they  are  off  by  90°.  For  small  R2, 
the  second  update  dominates.  As  a  result,  the  minor 
eigenvector  points  in  this  LOS  direction  and  the  major 
eigenvector  is  perpendicular  to  it. 

For  large  R2,  the  second  update  becomes  inconsequential 
and  the  first  update  dominates.  That  is  why  the  curves 
remain  close  to  ft  as  ft  varies  its  value. 

It  is  interesting  to  note  the  curves  for  R2  =  R,  seemly  have 
a  jump  of  90°  at  ft  =  120°  or  ft  -ft  =  90°.  This  is  actually 
due  to  the  switch  of  labeling  of  major  and  minor  for  the 
same  eigenvectors.  This  is  best  illustrated  in  Fig.  12  for 
eigenvalues  and  Fig.  13  for  eigenvectors  for  R2  =  Rj  as  a 
function  of  ft. 

As  shown  in  Figs.  10  and  11,  when  R2  and  Rj  are 
comparable,  the  minor  eigenvector  angle  lies  between  ft 
and  ft.  As  the  second  LOS  vector  ft  moves  away  from  the 
first  LOS  vector  ft,  the  minor  eigenvector  also  moves  away 
from  ft.  But  the  eigenvector  angular  displacement  lags  that 
of  ft  and  stays  somewhere  between  ft  and  ft. 

Based  on  the  simulation  results,  heuristic  curves  are  shown 
in  Fig.  14  for  analysis  and  design.  Fig  14(a)  shows  the 
approximate  eigenvalues  as  a  function  of  angular 
separation  between  the  two  LOS  vectors.  Fig  14(b)  shows 
the  approximate  eigenvector  angles  as  a  function  of  angular 
separation  between  the  two  LOS  vectors.  Linear 
interpolation  can  be  used  to  find  the  eigenvalues  and 
eigenvector  angles.  In  the  figures,  the  effects  of  Ri  and  R2 
are  ignored  but  could  be  added  with  different  curves  in  the 
same  plot. 

For  an  arbitrary  covariance  matrix  A  and  its  inverse  A"1,  the 
eigen-decomposition  is  written  as: 


A  =  UAUr  (24a) 

A-1  =  UA^lf  (24b) 

where 

U  =  [u,  •••  uj  (24c) 

A  =  diag([At  •••  2.J)  (24d) 

A”'  =  diag([^'  •••  A~'])  (24e) 


are  the  eigenvectors,  eigenvalues,  and  inverse  eigenvalues, 
respectively. 


The  matrix  A  and  its  inverse  A"1  can  be  further  written  as: 

A  =  ^T1uiuf  (25a) 

1=1 

A-1  =  ^Tr'iquf  (25b) 

!= 1 

Since  the  eigenvectors  are  orthonormal,  they  can  be 
considered  as  the  LOS  vectors  as  in  (22).  As  a  result,  the 
update  of  a  unity  circle  with  a  covariance  matrix  (or  a  full 
rank  observation  matrix)  can  be  readily  obtained  from  (23), 
leading  to  the  updated  ellipse,  which  is  re-oriented  by  the 
original  eigenvectors  (23b)  but  re-sized  by  the  updated 
eigenvalues  (23a).  This  thus  provides  a  theoretical 
justification  of  our  intuition  about  LOS  updates,  a 
geometrical  interpretation  of  combining  two  ellipses,  and  a 
computational  procedure  for  fusing  them. 

Since  the  trace  of  a  matrix  is  the  sum  of  its  eigenvalues,  we 
can  evaluate  the  geometric  dilution  of  precision  (GDOP) 
directly  from  eigenvalues.  For  the  case  of  updating  with 
two  LOS  vectors  where  ft  =  30°  and  Ri  =  1,  Figure  15 
shows  the  traces  as  a  function  of  ft  for  different  R2.  When 
ft  =  ft,  the  trace  is  2-f ,  which  is  the  largest  value.  The 
smallest  value  occurs  at  ft  -  ft  =  90°,  which  is  2  -  y  -  y2.  At 
other  ft,  the  trace  lies  in  between  following  a  sinusoidal 
curve.  Fig.  16  is  a  linear  approximation.  For  large  R,  the 
trace  does  not  change  much  with  ft  But  it  is  significant  for 
small  R.  As  shown  in  Fig.  15,  the  change  in  trace  over  A 9= 
90°  is  about  the  same  for  R2  =  0.1  to  10  with  R,  =  1.  This  is 
a  tradeoff  between  assigning  sensors  with  good  quality  vs. 
with  good  geometry  over  different  time  horizons. 

The  above  analysis  has  been  applied  to  the  general  case 
with  an  arbitrary  initial  error  covariance  [11].  Its 
relationship  with  GDOP  [13]  as  well  as  geometric  measure 
of  merits  (GMOM)  [4]  is  presented  in  the  context  of  sensor 
resource  management  for  layered  sensing  [12]. 

4  Conclusions 

In  this  paper,  the  use  of  range  and  bearing  measurements 
for  target  positioning  and  tracking  was  investigated  from 
the  geometric  point  of  view.  As  a  scaling  factor  to  the 
measurement  prediction  error  vector,  the  observation 
matrix  was  shown  to  be  made  of  line  of  sight  (LOS)  vectors 
from  ranging  sensors  to  target  and  unit  vectors 
perpendicular  to  the  LOS  scaled  by  range  from  bearing- 
only  sensors.  As  shown,  the  LOS  vectors  together  with  the 
sensor  measurement  error  covariance  determine  the 
eigenstmcture  of  position  estimation  error  covariance, 
shaping  the  error  ellipses.  This  geometric  insight  into  target 
positioning  and  tracking  with  range  and  bearing 
measurements  is  helpful  in  designing  an  efficient  resource 
management  strategy  for  distributed  sensors. 
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Fig.  8  Max  Eigen  over  A# for  Different  cr 


Fig.  9  Min  Eigen  over  A# for  Different  <j 
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Fig.  13  Max  and  Min  Eigenvectors  over  ts.0 


91 


Figure  14  Linear  Interpolation  for  Eigenvalues  and  Eigenvector  as  a  Function  of  A# 
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Fig.  15  Trace  vs.  A#for  Different  cr 
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